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Individual transmissions must be identifiable to distant 
receiving stations and must be synchronized so that no 
interference overlap occurs. A unique word sequence can 
serve both purposes. 

The basic .concepts of a satellite communications TDMA 
system are briefly reviewed, with emphasis on functional 
receiver requirements. A correlation detection scheme is 
hypothesized, from which false and miss-detection 
probability expressions are derived. criteria are then 
established for the selection of sequences suitable for use 
as unique words. Representative unique words of length 15 
through 20 are included, along with an explanation of the 
difficulties encountered in devising an efficient computer 
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I. INTRO DUC TION 



The possible use of time-division multiple-access (TDMA) 
in satellite communications has received a considerable 
amount of attention gin recent years. In a TDMA system, many 
stations use the same satellite repeater on a time shared 
basis. Transmitter timing must be synchronized so that only 
one station's signal accesses the repeater at any given 
time; and, further, so that successive signal bursts are 
closely spaced in order to make efficient use of the 
repeater capacity. Continuously changing ranges from ground 
stations to the satellite (Appendix B) make transmitter 
synchronization a non-trivial problem. Once the 
synchronization problem is solved, TDMA offers several 
advantages over current frequency-division multiple-access 
procedures [1,2]. Schmidt [3] outlined the proposed 
application of TDMA to the lliTEL5AT-4 series commercial 
communication satellites. 

A. TDMA SYSTEM DESCRIPTION 

As each station must transmit in timed bursts, 
transmissions are necessarily digital and provisions must be 
made for the sampling (and multiplexing) of the transmitted 
information channel (s). This process must be timed so that 
it is in consonance with the burst repetition freguency 
(BRF) and transmitted data rate. When analog-to-digital 
conversion is involved, the analog sampling rate must be 
egual to, or a multiple of, the BRF; a sample storage 
capability is required when the sampling rate is a multiple 
of the BRF. When information to be transmitted is in the 
form of a continuous bit stream, a buffer is required to 
store a given portion of the bit stream for transmission at 
the (higher) data transmission rate. In either case, the 
number of information channels that can be transmitted in a 
burst of given length is a direct function of the ratio of 
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transmission rate to information rate. Time-division 
multiple-access, then, is in effect a higher level form of 
time-division multiplexing. TDMA is a time sharing at RF, 
while TDM is a tj.roe sharing at baseband, before carrier 
modulation. The process of several transmitted signals time 
sharing a satellite repeater is conceptually similar tc that 
of several information channels time sharing the same RF 
carrier. 

Given that transmitted signals are neatly interleaved in 
their arrival at the satellite, they will be retransmitted 
in the same fashion and received on the ground without any 
interference overlap. A receiver must first recognize when 
a burst is of interest and then, in the case of multiplexed 
transmissions, select particular information channels within 
the burst and route them to appropriate end terminals. In 
order to do so, the receiver must be provided with certain 
control information. 

1 • Receiver F unctio nal Requ iremen ts 

Successive received bursts have been originated by 
different, in general widely separated, stations whose 
carrier frequencies are not coherent in frequency and phase. 
Therefore, prior tp the demodulation of each burst, the 
receiver local oscillator must be locked in frequency and 
phase to that of the incoming carrier using phase-locked 
loop techniques. (Use of differential phase-shift keying 
removes the carrier recovery requirement at the expense of 
error-rate performance £4].) Similarly, once carrier 
recovery has been accomplished, the receiver bit timing 
oscillator (clock) must be synchronized to the incoming bit 
stream. Bitzer £5] proposed a system design for 
establishing coherent bit timing among all stations and 
Suguri, Doi and Metzger [4] reported on the successful test 
of this type of system. In such a system, receiver and 
transmitter clocks are continuously adjusted to maintain 
coherency within acceptable limits, obviating the need for 
gross adjustments to the receiver clock preparatory to burst 
demodulation . 
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With the l,ocal oscillator and clock synchronized, 
demodulation of the incoming signal occurs. Embedded 
somewhere in the beginning of the burst must be information 
whiph the receiver can extract to obtain transmitting 
station identity, word (channel) synchronization and channel 
status. Demodula tion/demultiplexing of the information 
channels commences thereafter. 

2. Frame and Burs t Form at 

Satellite repeater access time is divided into 
frames, a frame being the length of time allotted for the 
transmission of one burst from every station. The frame 
length is typically 125 usee, or a multiple thereof, to 
conform to the Nyguist sampling rate for voice (8KHZ) [3,6], 
A station is assigned a time slot within the frame for the 
transmission of its burst, the length of which depends on 
the anticipated channel requirements. A guard time of 
typically 200 nsec or less separates the time slcts to 
preclude burst overlap. 

Each burst commences with a preamble, followed by 
the information channels. A typical format is depicted in 
Fig. 1-1. A fixed pattern of dummy bits is inserted at the 
beginning for the purpose of locking up the local oscillator 
and clock. A constant string such as 0Q00»©® for carrier 
recovery, followed by an alternating pattern such as 1010««« 
for bit timing recovery is typical. Following the recovery 
bits is a fixed sequence of bits unique to the transmitting 
station, the “unique word.'* More will be said about the 
unigue word below. The remainder of the preamble contains 
control signals and orderwire transmissions necessary for 
assignment and control of the information channels. 
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Figure 1-1 

Frame and Burst Format 



TDMA efficiency is a measure of the total time 
allotted to the transmission of actual information versus 
the total available repeater time. The preamble is wasted 
transmission time in this sense and considerable effort has 
been expended to reduce its length [4,5, 8, 9]. 

3 . Unigue W or cl Func tions 

To the extent that it is recognizable tc the 
receiver, the unigue word can be used to identify the 
originating station. It may also be used as a timing 
reference to establish word synchronization, resetting a bit 
counter which properly routes the control and orderwire 
portions cf the preamble and indicates the beginning cf the 
first information channel. 

It has been assumed up to this point that all 
transmitting stations are synchronized so that no burst 
overlap occurs. Nothing has been said about how 
synchronization is acquired or maintained. The usual method 
of doing sc is to use the unique words as timing marks 
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within the frame. One station is designated as the master 
station and its unique word used as the fundamental timing 
reference for the frame. All other stations compare the 
time difference between their own unigue word and the 
reference, as received from the satellite, and adjust 
transmitter timing as necessary, one bit per frame, to 
maintain proper separation. This problem of transmitter 
synchronization is of interest herein only in that the 
unigue word plays a vital roie, and will not be further 
pursued. references 3,7,9 and 10 offer a. detailed 
description cf the various factors which must be considered 
and provide insight as to how a synchronization package 
might be instrumented. 

B. DETECTION AND SELECTION OF UNIQUE WORDS 

A TDMA system of the type outlined above can be 
operated successfully only if the receiver is able to detect 
the receipt of a unigue word and differentiate one from the 
other. Some type of correlation detector is inferred. 
Proposed detector design features and operation are 
presented in Section II. Then, in Section III, criteria are 
established for the selection of unigue words which yield 
acceptable miss-*detection and false detection probabilities. 

Interest in the subject was motivated primarily by the 
discussions of Schrempp and Sekimoto [11]. Much of the 
detector design philosophy is borrowed from their work, and 
their approach to unique word selection has been 
instrumental in obtaining an appreciation for the 
fundamental requirements which must be considered. The 
selection criteria set forth in Section III and the 
subsequent computer search results represent the major 
products cf this research effort. 
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II. DET ECTI ON OF UNIQUE WORDS 

A unigue word is a fixed binary sequence of length N. 
Given a knowledge of the elements of the sequence, detection 
can be accomplished .by comparing received bits, N at a time, 
with a replica of the seguence stored in the receiver; i.e., 
by using correlation techniques. The probability of missing 
the unique word (miss-detection) and the probability of 
receiving the sequence at random (false detection) are 
fundamental considerations in the evaluation of detector 
ef f ectiveness. 

A. EASIC DETECTOR DESIGN 

A digital correlation detector of the type depicted in 
Fig. 2-1 is postulated. The element values of the unique 
word are set semj.-permanently in a register and the 
demodulated tit stream is clocked through the shift 
register. The corresponding elements of the two registers 
are added, mcdulo-2, once each bit period and the results 
are summed algebraically for input to the threshold network. 

The mcdulo-2 operation is defined by, 

0 © 0 = 0 
0 9 1 = 1 
19 0=1 
19 1=0 

And the input to the threshold network is, 

R = 2 H (s « U ) (2.1) 

r= 1 i i 
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Figure 2-1 

Digital Correlation Detector 



R therefore equals the number of disagreements between S 
and U. The threshold network is designed so that a 
detection pulse is produced whenever there are h or fewer 
disagreements; that is, whenever R<h. The setting of a 
threshold permits the detection of unique words which 
contain up to h transmission errors and thereby reduces the 
miss-detection probability. The probability of false 
detection increases in doing so, however, due to the 
increased opportunities for a random match-up. 

It should be noted that R, as mathmatically defined, is 
in fact the inverse correlation function of S and U. When 
the elements of S and U all match, R=0; and when there is 
complete disagreement, R=N. Alternative detector designs 
can be envisioned which would produce the true correlation 
function at R; e.g., construct adders that produce the 
complement of the modulo-2 sum and design the threshold 
network so that a pulse is produced whenever R>N-h. Of 
course, in the end, the voltage levels associated with 'O' 
and '1' are what determine the actual "value" of R. 
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Whichever of the various alternatives are chosen, the end 
result is the same: a detection pulse is produced whenever 
there are h or fewer disagreements between S and U. 

The capability of receiving more than one transmitting 
station is provided by a parallel bank of detectors, each 
set to the unique wprd of one of the transmitting stations. 
The received bit stream is passed simultaneously through 
each detector and the resulting detection pulses are used to 
control receiver functioning, as described in Section I. A. 3. 

B. STEADY-STATE APERTURE GATING 

The probability of false detection is reduced 
significantly if the detector can be gated so that it begins 
looking for a unique word just prior to its receipt, and is 
otherwise inactive. The detection pulse itself may be used 
to generate such a gating signal under steady-state 
conditions; that is, once the unique word has been initially 
apguired . 

A unique word is transmitted every I seconds, once each 
frame, except when a correction is required to maintain 
proper position within the frame. If corrections are made 
at the rate of one bit per frame [9,10], the elapsed time 
between unigue word ftransmissions is bounded by T + b, where b 
is the period of one bit. The time separation of successive 
received unigue words is further affected by relative motion 
between the satellite and the respective transmitting and 
receiving stations. Let d, in bits per second, be the 
maximum expected doppler induced drift of successive unigue 
word detection pulses. The maximum doppler shift in one 
frame period is then ±Tbd and the elapsed time between 
successive receipts of a unigue word is bounded by, 

T = T ± (b + Tbd) 
r 

T = T ± b ( 1 + Td) (2.2) 
r 
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This relationship can be expressed in terms of bit 
periods by dividing through by b. Doing so establishes a 
limit on the number of bit periods separating successive 
receipts of a unigue word. 

N = T /b = T/b ± (1 + Td) 
r r 

N = n ± (1 + Td) 
r 

where n is the number of bits in one frame. 

Defining D as tie maximum displacement (in bits) cf the 
unigue word from its "periodic" position, 

D = 1 + Td (2.3) 

N = n ± D (2.4) 

r 

Upon receipt of a unigue word, the detector is turned 
off, gated back on T-bD seconds later, gated off by the next 
unigue wolo, etc. False detections can therefore result 
only from a random match-up within the aperture immediately 
preceding the unigue word. 

The guestion remains as to what to do when the unigue 

word is missed due to errors occurring in transmission. 

Once the detector is gated on, is it left on until a 

detection pulse is received (in which case a missed unigue 

word causes an immediate transistion to the acguisition 

mode) ; or, is the detector gated off at time T+bD and then 

1 

gated back on at time 2T-bD, one frame later? In the 
latter case, how many successive misses are reguired before 
cancelling the aperture and reverting to an acguisition 
mode? The answer depends for the most part on particular 



1 

It is assumed that the aperture width remains constant 
and would not be expanded as successive misses occur. 



17 



system specif icatioas. When a unigue word is missed, all 
information contained in the burst is lost. Having missed 
it once, that is having received the unigue word with more 
than h errors, there is a given probability that it will 
contain more than h errors the next time. Calculation of 
the probability of a unigue word containing more than h 
errors two or more times in succession is not enough 
however. The doppler induced burst drift rate must also be 
considered. 

if the drift rate is such that the unigue word is 
received outside the aperture, it will be missed, even if 
error free. On the pther hand, reversion to an acguisition 
mode may involve a considerable penalty in lost information 
due to a long average acguisition time. (See Section II. C 
below.) Thus, a desire to avoid the acguisition mode must 
be balanced against the diminishing probability (as 
successive misses occur) of receiving the unigue word in the 
aperture mode. As stated earlier, this is a systems 
dependent decision, and one that can best be reached 
empirically in an operational environment. For present 
purposes, the rather safe assumption will be made that at 
least two successive misses are reguired before a transition 
from the aperture to acguisition mode is effected. A 
detection pulse turnis off the detector and sets a counter 
which gates the detector on T-bD seconds later. If no 
detection pulse occurs, the counter turns off the detector 
at time T+bD. The aperture is of width 2D+1 (bits) . 

The maximum burst drift which may be expected when using 
a geosynchronous satellite is computed in Appendix B. 
Eguation B.5 can be substituted into eguation 2.3 to obtain 
values for the maximum displacement of a unigue word over 
one frame period. Table 2-1 contains representative 
calculated results for the shipboard platform discussed in 
Appendix B as well as for a platform with a relative 
velocity on the order of 1000 knots. 
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Table 2-1 

Bit Rate Doppler (d) and Maximum Displacement (D) 
of a Unique word Over One Frame Period for Given 
Values of Station/Satellite Relative Motion (v) , 
Bit Frequency (f) and Frame Period (T) 



C. ACQUISITION GATING 

The steady-state operation of a unique word detector was 
described in the preceding section. Here the process of 
acquisition is considered. When a receiver is activated, 
each detector must acquire the associated unigue word within 
the TDHA frame. Likewise, a detector must be switched to' 
the acquisition mode whenever steady-state synchronisation 
is lest. Acquisition, or re-acquisition, may be effected by 
simply passing the received bit stream through the detector 
until a detection pulse is received. Then, as in the 
steady-state mode, tie detection pulse switches the detector 
off and starts an aperture gate counter which in turn 
switches the detector back on in the next frame. If a 
second detection pulse is received within the aperture, the 
detector is transferred to the steady-state mode. 
Otherwise, the acquisition process continues. 

When in the acquisition mode, the bursts of all active 
stations are being passed continuously through the detector 
and the probability of false detection is thus relatively 
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high. Should a false detection occur, a false aperture is 
set and the next actual unique word detection pulse will be 
missed.. Even when the detection is valid, there is a 
non-zero probability that the unigue word will be missed in 
the next frame. In either case, acquisition, and transfer 
to the steady-state mode, is delayed for at least two frame 
periods. The actual acquisition time may thus span several 
frame periods. 

It is possible to construct a detection pulse storage 
network which would look at detection pulse locations over 
two or more frame periods and use majority logic to set an 
aperture corresponding to the true detection pulse. . 
Although acquisition time would be improved, implementation 
could prove relatively costly because each detector would 
require its own network. An alternative is proposed, which 
makes use of the same aperture network employed in the 
steady-state mode to gate the detector and reduce the 
probability of false detection. 

1 • frame Bef erence Gating 

Each burst has an assigned position within the TDMA 
frame. It is known that a particular burst will be received 
a given number of bits after the receipt of the reference 
unigue work, plus or minus the guard band cushion. ' The 
guard band is set so that it prevents burst overlap; that 
is, it represents the maximum expected displacement of a 
burst (or unigue word) from its nominal position. It is 
therefore proposed that, during acquisition, the 
aperture/counter network of a detector be keyed by the 
receipt of the reference unigue word. The aperture is set 
so that the detector is gated on at the assigned position of 
the corresponding unique word, minus the number of bits in 
the guard band, B. If the unigue word is missed, the 
detector is switched off B bits after the nominal unigue 
word position. (The aperture width is 2B+1.) If the unique 
word is detected, the aperture is set for the next frame. 
Two consecutive detections key a transfer to the 
steady-state mode. 
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In a demand access system, where at least some of the 
stations are not assigned a permanent time slot, use of 
frame referecce gating is feasible only if all stations can 
be provided with time slot activation and deactivation 
information, or can otherwise predict where a given 
station's burst may appear within the frame. Whether cr not 
it is possible to provide advance activation notification to 
all stations depends on how a part time station gains access 
to a time slot. If access is centrally controlled 
(a utomafcically) , then assignment signals from the 
controlling station can be copied by all stations and used 
to set the acquisition gate. If a station gains access by 
seizing a vacant time slot, a receiving station must have 
other means of predicting where the burst will occur. The 
results of the computer search for unique words, reported in 
Section IV, suggest a procedure for meeting this 
requirement. 

For a given unique word length, several different sets 
are available which have good cross-correlation properties. 
Unigue words from any one set could be assigned to certain 
stations which collectively have exclusive use of a block of 
consecutive time slots within the frame. As many such 
blocks as necessary could be made available to different 
groups of stations.. The number of time slots assigned to 
each group would be based on collective useage requirements, 
and could be varied on a scheduled or demand basis. (Demand 
access implies that ±he number of time slots would be less 
than the number of stations.) Those receiver detectors 
corresponding to the demand access transmitters could then, 
when in the acquisition mode, be gated on at the beginning 
of each burst uithi-n the demand access sub-frame. The 
probability of false detection will be somewhat higher than 
for the single gate procedure proposed above; but, with low 
cross-correlation between unigue words of a set, will be 
much lower than for ."wide open" acquisition. 
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III. UNIQ UE WORD SELECTION CRITE RIA 

Having outlined the basic detector design concepts, a 
basis has been formed for the selection of sets of unique 
words. A detailed analysis of the attendant detection 
probabilities is presented in Section C. Simply stated, the 
goal is to find uni.gue words which yield low miss-detection 
and false detection pr obabili ties. Criteria will be 
established in this section for doing so. 

A. LENGTH 

The fewer bits in a unique word, the lower the 
miss-detection probability. Conversely, the longer the 
unigue word, the lower the false detection probability. The 
use of a detector threshold and aperture nullifies to some 
extent the importance of these opposing considerations, but 
in a gross sense the.y can never be eliminated as factors 
influencing the ultimate choice of unique words. In fact, 
length and the threshold setting are the cnly parameters 
available for altering the miss-detection probability. 

From the standpoint of efficiency, it is desirable to 
keep unigue words as short as possible. The longer each 
station's unigue word, the less time is available for the 
transmission of information. 

E. AUTOCORRELATION CRITERIA 

Operation of the detector in the steady-^sta te mode is 
described in Section II. B. Loss of channel information will 
result if (1) the unigue word is missed, or (2) a false 
detection occurs in the aperture prior to the unique word. 
Either may lead eventually to a drop-sync condition. As (1) 
and (2) are equally damaging, it would seem appropriate to 
select unigue words in such a way that their combined effect 
is minimized. Under the conditions described in Section 
II. E, it is largely a matter of u ng ua ntif iaLl e chance 
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whether a miss-detection or a false detection (within the 
aperture) is the least desirable occurrence. Accordingly, 
the following condition was established: if the received 

unigue word contains h or fewer errors, the false detection 
probability must be zero within the aperture. A false 
detection pulse will not be produced unless the unigue word 
is going to be missed anyway. 

As it is possible for a unigue word detection pulse to 
occur on the last cycle within the aperture, and as the 
aperture is of width 2D-M, tho foregoing condition can be 
met only if the probability of false detection is zero for 

the 2D cycles preceding the unigue word. This condition can 

in turn be satisfied only if the autocorrelation function of 
the unigue word is such that, 

Rj(t) < N - 2h - t ; 1<t<2D (3.1) 

where t is an integral time variable, the basic unit of 

which is one bit period, and, 

N-t 

B(t) = ? U(i)©0(i + t) 

1=1 

This condition is presented graphically in Fig. 3-1. 
R (t) is egual to the number of agreements between the 
(error-free) N-t displaced unigue word bits and the detector 
setting. Up to h errors are permitted in the unigue word 
under the prescribed conditions and it is possible that 
these errors may b>e such as to produce additional 
agreements, giving a total of R(t)+h agreements from the 
displaced unigue word alone. There are t non-unigue word 
bits in the detector. If all of them match with the 
detector setting, a combined total of R(t)+h+t agreements 
results. As N-h or more agreements are reguired to produce 
a false detection, the probability of such occurrence is 
zero if R (t) +h*t<N-h, and the validity of ineguality 3.1 is 
established. 
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FIGURE 3-1 

AUTOCORRELATION CRITERIA 



In the following section, a condition for the 
cross-correlation function will be established. This same 
condition must be met by the autocorrelation function 
outside the aperture in order to meet specified limits on 
the probability of false detection during acquisition. It 
is listed here, in advance, in order to complete the 
autocorrelation specifications. 

R (t) < N - h - t ; 2D<t<N-h (3.1a) 

R (t) = 0 ; N-h< t <N- 1 (3.1b) 

N 

There are 2 sequences of length N from which to choose 
unigue words. However, symmetry considerations can be 
invoked to reduce the number which must actually be checked. 
First, because the m,od-2 summing operation is commutative, a 
seguence and its complement, obtained by inverting the value 
of each bit, have the same autocorrelation function. 
Second, because autocorrelation is an even function, a 
sequence and its transpose (or mirror) , obtained by 
inverting the order of the bits in the sequence, have the 
same autocorrelation function. Therefore, once a seguence 
is checked it is unnecessary to check its complement, 
transpose or complement of its transpose. The number of 
required checks is npt reduced by a factor of four, as it 
would appear at first glance, however, because the transpose 
and complement transpose pair may be identical tc the 
original seguence and its complement. The required number 
of checks can be determined by consideration of the effect 
of the complement and transpose operations. Although the 
results were not used directly in the computer search 
routine, they do provide an indication of search efficiency, 
and are presented here as a matter of peripheral interest. 

In order for the transpose or complement transpose to be 
identical to the original seguence, the second [N/2] bits of 
the seguence must be the transpose or complement transpose 
of the first [ H/2 ] bits, where [ N/2 ] equals the largest 
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integer iess than or egual to N/2. (If N is odd, the middle 
bit of the seguence does not change position under the 
transpose operation,.) The first [N/2] bits may be arranged 
[ N/2] 

in 2 different ways, each arrangement having an 

associated transpose and complement transpose in the second 

[N/2] . . N 

half. Therefore, 2x2 of the original 2 sequences have 

redundant (identical to original seguence) transposes or 

complement transposes, and, for this subset, the required 

number of checks can be reduced by a factor of two only (by 

pairing a sequence with its complement),. There are 
N [ N/2 ] + 1 

2 -2 sequences remaining to be considered. These can 

be grouped by fours (sequence, complement, transpose, 
complement transpose) to reduce the number of reguired 
checks by four. The total number of sequences of length N 
which must be checked in search of unique words meeting the 
autocorrelation condition is therefore. 



Acquired Checks 



N 

2 



2 [N/2]+1 

4 



2 1 "/ 23 



N-2 [ N/2 ]- 1 

= 2 +2 

The computer program used to accomplish the 

autocorrelation sort is listed in Appendix D. 

C. CROSS-CO ABEL ATI ON CRITERIA 

The c^oss-correlation between the unique words used by 
stations sharing the TDM A frame is a crucial factor 

affecting the false dietectic probability during (ungated) 
acquisition. A cross-correlation characteristic that yields 
a high probability of false detection is particularly 
unsatisfactory because of the periodicity of its occurrence. 
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In such cases, the detection aperture may be falsely set 
time and time again pver several consecutive frames. 

The probability expressions derived in Appendix C were 
used as a basis for the establishment of the following 
cross-correlation cpndition; the probability of false 
detection must be zero unless the received unique word 
contains cne or more errors. This condition also applies to 
the head and tail of the autocorrelation function, outside 
the aperture. The cross-correlation condition is thus 
similar to the autocorrelation condition, differing only in 
that the unigue word must be error free, vice h or fewer 
errors, to guarantee zero false detection probability. 
Using B ( t) again to xepresent the correlation function, the 
condition cn R (t) is (analogous to inequality (3.1)) 
R (t) <N-h- J 1 1 . As .the correlation function is always 
non-negative, this condition can be met only when |t{<N-h. 
For higher values of |t| the condition must be relaxed to 
R (t) =0. Summarizing, sets of possible unique words, having 
passed the autocorrelation test, are retained only if the 
cross-correlation be.tween all words of the set satisfy, 

R(t) < N - h - 1 1 1 ; 1 1 1 <N-h (3.2a) 

R (t) = 0 ; N-h< 1 1 |<H-1 (3.2b) 

The conditions placed on R (t) can be used to further 
refine the probability expressions of Appendix C. It will 
be assumed that bits of fixed value (at least N-1 of them) 
precede the unique word and that bits of random value 
follow. This assumption is based on the system 
specifications outlined in Section I, and, although not an 
absolute necessity at this point, is made in order to reduce 
the number of situations which must be considered. Other 
situations, when and if occurring, may be handled in a 
completely analogous fashion. 

First, the probability of false detection in the fixed 
bit region preceding a unique word (a unique word other than 
the one for which the detector is set) will be considered. 
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FIGURE 3-2 

CROSS-CORRELATION CRITERIA 



The point where the unique word is completely contained in 
the detector defines t=0. At t=-(N-1) the first bit cf the 
unigue word enters the detector. Eguation (s) C.8 apply for 
-(N-1)<t<0, with m=N. As in Appendix C "number of 
agreements" refers to matches between the detector setting 
and the value of the fixed bits as transmitted. The cumber 
of agreements, k in eguation C.8, is egual to R (t) , the 
number of matches within the unigue word bits, plus the 
number of matches within the remaining Jt| bits. If all of 
the latter match. 



k = R(t) <- |t | 

Therefore, the cross-correlation condition 3.2a 
(R (t) <N-h- 1 1 1 ) is equivalent to the eguation C.8 condition, 

k + h < N 

The cross-correlation condition 3.2b in itself does not 
guarantee that k + h<N. It is necessary to impose the 
additional systems dependent condition that fewer than N-h 
agreements can be permitted between the non-unigue word bits 
for - (N- 1)<t<- (N-h) . Doing so, the probability of false 
detection at any given step is (from eguation C.8e), 

P{A(t)} = 

N-h-k h+k 

C(N-k,N-h*Jc) p (1~p) .* - (N- 1 ) <t<0 (3.3) 

If k takes on its maximum value (N-h) , the probability 
of false detection becomes, 

N-1 

P{A(t)} = (h+1) p (1-p) ; - (N-1) <t<0 (3.4) 

max 

The probability of false detection in the random bit 
region following the unique word is also given by 
eguation (s) C.8. In this case m equals the number of unique 
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word bits in the detector (there are no other fixed bits) 
and the number of fitted matches, k, equals R (t) . 



m = N - |t | 
k = E (t) 

Therefore, the cross-correla tion condition 3.2a 
(R (t) <N-h- | t | ) is equivalent to the equation C.8 condition, 

k + h < m 

And the probability of false detection at any given step is 
given by equation C.Bc (with m and k given by above) . 

P{A(t)} = 

N-m m-k-h k+h 

(1/2) C (m-k, m-k-h) p ( 1 — p) ; 0<t<H-h (3.5) 

When k takes on its maximum value (R (t) =N-h-t- 1) , 
equation 3.5 becomes, 

t N-t-1 

P{A(t)} = (1/2) (k+ 1) p (1 — p) ; 0<t<N-li (3.6) 

max 

For N-h<t<N- 1 , m=h,h-1 ,h-*2,»®«,1 and k=R(t)=0, giving 

m = N - t 

and , 

k + h > m 

The probability ,of false detection is (from equation 
C. 8b) , 



N-m ra h-m 

P{A(t)} = (1/2) (1-p) 2 n c (N-m, j) (3.7) 

3 = 0 

t N-t h-N + t 

= (1/2) (1-P) .2 C (t , j) ; N-h<t<K-1 

3 = 0 
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Equation 3.7, unlike equations 3.3 and 3.5, is 
independent of parameter variation [i.e. independent of k, 
fi (t) ] and thus "maxijnum P{A(t)}" has no meaning. 

Equations 3.3, 3.5 and 3.7 may be used to calculate the 
probability of false detection as each unique word in the 
frame is passing through the detector. Equations 3.4, 3.6 
and 3.7 provide a worst case limit for these probabilities 
(when unique words are chosen using the cross-correlation 
criteria 3.2a and 3.2b). These individual fixed-bit event 
probabilities are the ones needed, among others, in equation 
C.10 to obtain an upper bound for the probability of one or 
more false detections (F) over the period of one frame. Or, 
over a narrower time span, they can be used to compute an 
upper bound on the probability of one or more false 
detections during the 2N-1 detection events surrounding a 
given unigue word; that is, to compute. 



F ! 



* 



N-1 



t=-N+ 1 



f * / A. \ > 

* t* V W i 



F* provides a means of obtaining a relative measure of 
the goodness of different unigue word choices without having 
to resort to the more lengthy calculations required by 
equation C.10. Substituting from equations 3.4, 3.6 and 

3.7, 

N-1 N-h t N-t-1 

F' < N[ (h+1) p (1-p) ] + | (1/2) (h + 1) p ( 1 P ) 

N-1 t N-t h+t-N 

+ 2 , (V2) (1-P) 7 c (t, j) (3.8) 

t=N— h- 1 3=0 

The computer program used to check cross-correlation is 
listed in Appendix D. Each sequence which passes the 

autocorrelation test is cross-correlated with every other 
such sequence, yielding pairs of acceptable unique words. 
In order to form maximal length sets of acceptable unique 
words, it is necessary to perform an additional matching 
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routine. For example, if UW 1 cross-correlates well 
(matches) with UW2* UW3 and UW4, then a set of three is 
obtained if UW2 also matches with UW3 or UW4, or UW3 matches 
with UW4; a set -of four is obtained if UW2 matches with both 
UH3 and UW4 a nd UW3 matches with 0W4. Various matching 
algorithms were programmed, none of which performed in a 
completely acceptable fashion. At N=17, the number of 
checks which must be performed becomes so great that, for 
the algorithms used, computation time runs into days. In 
order to obtain results for N=17 and higher, various 
expedients were employed, as described in the following 
section, to reduce tie required computer time. 
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IV. PM3MIATI0N OF RESULTS 



The SES-9300 digital computer was used to perform a 
correlation sorting routine patterned after the criteria 
established in the preceding section. At the outset, it was 
realized that a significant amount of machine time would be 
reguited to obtain results; but, even so, the magnitude of 
the problem was not fully appreciated until program 
execution was initiated. The SDS-9300 was chosen over the 
IBM-360 because of the local availability of multi-hour 
blocks of computer time, realizing that its limited on-line 
memory (32K) might prove restrictive. 

The choice of detector parameters and the essential 
findings of the computer search effort are covered in the 
following paragraphs.. Selected sets of unigue words are 
listed in Appendix A. 

A. CHOICE OF DETECTOR APERTURE AND THRESHOLD 

The motivation for use of a steady-state aperture 
technigue is provided in Section II. B. As the aperture 
width must be at least 2D+1, where D is the maximum 
displacement of a unigue word in one frame period. Table 2-1 
leads one to a choice of a four or five bit width. An odd 
width is preferred in order that the aperture counter can be 
run directly off the receiver clock; i.e., so that detector 
activaticn/deactivation occurs at the beginning rather than 
in the middle of a bit cycle. An aperture width cf five 
bits is thus the mimimum acceptable, and yet at the same 
time provides sufficient margin to cover unigue word drift 
over more than one frame in the event of a miss-detection . 

The choice of a detector threshold (h) is determined by 
the setting cf an acceptable limit for the miss- det ection 
probability (Q) , which in turn is a function of the bit 
error rate (p} . For this purpose, p was assigned a value of 
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-a 

10 [ 1 1 J as being d;he worst-case bit error rate that could 

be tolerated on a satellite communication channel. 

Eguations C.1, C.2 and C.3 were used to compute Q and the 

mean time between miss-detections (E) for various values of 

N and h. The resalts are listed in Table 4-1. The times 

given for H assume a frame period of 125 usee. 

Acceptable values for Ti occur for h=2 and higher. The 

fact that the probability of false detection increases with 
h dictates the selection of the minimum acceptable 

threshold; viz. h=2. 



N 


h 


0 


2 


15 


0 


1.50x10 


83 


E3GC 




1 


-6 

1.05x10 


2 


min 




2 


-10 

4.55x1 C 


76 


hr s 




.3 


-13 

1.36x10 


29 


yrs 


17 


0 


-3 

1.70x10 


74 


msec 




1 


-6 

1.36x10 


92 


sec 




2 


-10 

6.79x10 


51 


hrs 




3 


-13 

2.38x10 


17 


yrs 


20 


0 


-3 

2. 00x10 


62 


msec 




,1 


-6 

1 . 90x10 


66 


sec 




2 


-9 

1 . 14x10 


30 


hrs 




3 


-13 

4.84x10 


8 


yrs 



Table 4-1 

Miss-Detection Probabilities 
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B. UNIQUE WORD SYMMETRIES 

Certain characteristics of the sequences which pass the 
autocorrelation test were used to simplify the 
cross-correlation sort. These same characteristics also 
provide a degree of flexibility in applying the results to 
particular operating systems. 

A seguence can pass the autocorrelation test only if its 
autocorrelation is .zero for the first and last three values 
at the head and tail of the distribution (with h=2) . 
Therefore, the first three bits of the seguence must be all 
zeros or all ones, and the last three bits must all be of 
opposite polarity. The need to check less than one-’half of 
the total possible sequences was explained in Section III.B. 
As the autocorrelation sort looks at sequences in ascending 
order of numerical value, sequences passing the test 
necessarily begin with three zeros. (The search is 
exhausted prior to reaching sequences which have a one as 
the first bit on the left.) When a potential unique word 
(UW) seg uence passes the autocorrelation test it is known 
that its complement (CUN) , transpose (TUN) and complement 
transpose (CTUW) will also pass the test. The latter were 
computed and stored along with the UW. For example, the 
following sequences ,of length 14 pass the autocorrelation 
test: 

UW CUW TUW CTUW 

00000101100111 11111010011000 11100110100000 00011001011111 
00001011001111 11110.100110000 11110011010000 000011001C1111 

The cross-correlation condition is identical to that of 
autocorrelation on the head and tail of the distribution. 
That is, the cross-correlation between sequences must be 
egual to zerc for the first and last three values. It is 
obvious from the example above that cross-correlations 
between the UW/CTUW sequences and the CUW/TUK sequences do 
not meet this condition. This is true in general. Further, 
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because or the symmetry inherent in the complementing 
operation, (1) if UW1 cross-correlates well (matches) with 
CTUW2, then CUW 1 matches with TUW2 and (2) if UW1 matches 
with UW2, then CTUW1 matches with CTUW2, etc. Therefore, in 
checking cross-correlation it is only necessary to check 
each UW against other UW’s and CTUW’s. Each set of unique 
words obtained by this process can be converted, by 
inspection, to a second set by taking the complement of each 
unique word. This provides a dual flexibility in applying a 
set to a particular operating system with its own peculiar 
arrangement of fixed value bits preceding the unique word. 

C. UNIQUE WCBD SETS 

As explained above, the sequences from the classes 
designated CUW and TUW need not be, and were not, considered 
in the cross-correlation search for sets of unique words. 
Each set extracted from the UW/CTUW seguences has a 
companion within the CUW/TUW seguences, obtainable by 
inspection. In order to accurately reflect the 
computational requirements , the reported results which 
follow thus are drawn only from the former. That is, such 
quantities as "number of seguences passing the 
autocorrelation test." and "number of sets of length X" refer 
to actual computer results. The numbers can be doubled by 
including the complemented sequences. Computer program 
descriptions are included in appendix D. 

Two results of primary interest were obtained for each 
value of N. First, the unique word which cross-correlated 
well with the maximum number of seguences was recorded, 
along with the corresponding matching seguences. These are 
referred to as "gated sets." These seguences may be used as 
unigue words in a gated acquisition system where it is 
necessary that only the reference unique word 
cross-correlate well with the other unique words. That is, 
as all detectors except the one set for the reference are 
gated, only the latter must contend with non-associated 
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unigue words during acquisition, and is thus the only one 
for which cross-correlation properties are of concern. 
Second, sets of unique words were obtained which have the 
property that every word in the set cross-correlates well 
with every other .word in the set. Sets of this type are 
necessary for ungated acquisition, and are referred to 
hereinafter simply as n sets." 

With h=2 and a five bit aperture width, no sequences 
passed the autocorrelation test for N=13 and below. At 
N-14, four sequences passed the autocorrelation test, but 
none of them could be paired in cross-correlation. At N=15, 
36 sequences passed the autocorrelation test. A gated set 
of length 13 and a total of 16 sets of length three were 
obtained. The total computer search time was 25 secs. 

At N=16, 150 sequences passed the autocorrelation test, 
and the largest gated set was of length 66. A total cf 56 
sets of length six were found. The computer search time was 
3 min, 50 secs. 

The difference in computational requirements between 
N=15 and N=1o was the first indication of the magnitude of 
the time problem which would be encountered in dealing with 
higher values of N. At. N=17, the number of 
cross-'-correlation pairs increased to the point where it was 
also necessary to re.vise the final matching algorithm to 
prevent memory overflow, compounding the time problem. Kith 
444 sequences passing the autocorrelation test, a gated set 
of length 272 was obtained. Several sets of length nine 
were found using 4+ hours of computer time, at which point 
program execution was terminated. It is estimated that 
approximately 65 hours would be reguired for a complete 
search with the algorithm being used. 

It was obvious that it would be futile to continue the 
search tc higher values of N using the same procedure. 
Attempts at developing a more efficient algorithm were 
unsuccessful. The 32K on-line memory of the SDS-9300 was a 
serious constraint, hut even without it, the anticipated 
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execution time for the algorithms envisioned far exceeded 
the computer time that could be made available. A simple 
subterfuge was eventually employed to obtain partial results 
for higher values of N. 

The number of sets of a given maximal length which were 
being found for higher values of N vere far greater than 
needed in order to offer an acceptable degree of flexibility 
in selecting a particular set for use in an operational 
system. Many of the sets differed only in the replacement 
of cne or two unique words, and, because of the nature of 
the searching process, tended to be uncovered within the 
same time frame. alien the program was terminated prior to 
completion cf a search, as at N=17, not only was it unknown 
whether or not sets pf greater length than found actually 
existed, the oqes that vere found were for the most part 
different combinations of the same slightly larger set of 
unigue words. Given the fact that a partial search must 
suffice, it seemed appropriate to effect a reduction in the 
number of seguences which passed the autocorrelation test to 
the point that the cross^correlation search could be run to 
completion. Doing so again creates uncertainty as to 
whether cr not sets of maximal length have been found. 
However, given the excessive number of sets being obtained, 
there is some hope that some of them might still appear when 
the autoccrrelation basis is reduced. 

Increasing the effective detector aperture creates a 
more strict autocorrelation condition and thus reduces the 
number of seguences passing the autocorrelation test. A 
mathematically artificial increase affects detector 
performance cnly in that the probability of false detection 
is reduced (unnecessarily) in the region immediately outside 
the real aperture. The effect of a varying aperture was 
first tested using N=15 and 11=16, for which complete results 
had been previously obtained. For N=15, the maximal length 
set remained at three for aperture widths of three through 
six bits. The number of such sets totaled 121 for a three 
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bit aperture versus six for a six bit aperture. For N=16, 
the maximal length remained at six for aperture widths of 
three through six bits, with corresponding totals of 724, 
345, 56 and 8. Four sets of length five were obtained for 
an aperture bidth of seven bits. 

Indeed, up to a point, the maximal length is independent 
of aperture variation. Partial results were obtained for 
values cf N equalling 17, 18, 19 and 20 using this 
technique. The aperture was set high and then reduced in 
increments until search time became excessive. The results 
so obtained again fell short of expectations in that much 
computer time was st^ll required to complete the partial 
search. For a given number of sequences passing the 
autocorrelation test, it was expected that search time would 
be of the same order of magnitude for the different values 
of N. In fact, search time increased rapidly with N. For 
example, with a five bit aperture, 150 sequences passed the 
autocorrelation test at N=16 and the total search time was 3 
min, 50 sec. With a 16 bit aperture, 142 sequences passed 
the autocorrelation test at N=20 and the search was only 
approximately 1/30th complete when terminated after two 
hours. The imposition of a stricter autocorrelation 
criteria produces a basis of the same order as for a lower 
value of N, but the cross-correlation properties of these 
sequences is such that many more combinations of match sets 
exist. That is, any given unique word cross-correlates well 
with many more sequences in the N=20 basis than in the N=15 
basis. 

Data on the numbers of sets found with h=2 for various 
values of N and various aperture widths is presented in 
Table 4-2; "W" is the aperture width and ''NAUTO" is the 
number of sequences passing the autocorrelation test. An 
asterisk denotes partial results obtained from an incomplete 
search. Table 4-3 contains similar data for h=1 and h=3, 
and serves to illustrate the effect of a varying threshold 
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on the numbers involved, 
maximal length sets .(for h=2) 



A representative listing 
is contained in Appendix A. 



of 
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N 
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6 
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95 
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6 
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11 
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10 
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102 


10* 


unkn 




9 
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19 


15 


36 


26 


5 
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14 
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83 


9* 
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13 
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20 


16 
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15 
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306 
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Table 4-2 

Sets Obtained For h=2 
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3 


19 
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0 


0 
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20 


7 


14 


0 


0 
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5 


,76 


17 


4 
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Table 4-3 

Sets Obtained For h=1 and h'=3 



D. FALSE DETECTION PROBABILITIES 

The probability of miss-detection is, for constant h, 
the same for all seguences of a given length. 
Representative number's are contained in Table 4-1. Given an 
acceptable combination of length and detector threshold, it 
is the relative values of false detection probability that 
then determine whether or not a particular sequence, or set 
of seguences, is suitable for use as a unique word. The 
auto- and cross-correlation criteria were established with 
this in mind. As a point of final analysis, the equations 
of Appendix C were used to calculate the false detection 
probabilities. Representative results are included in the 
following paragraphs. 

1 . Probability of Fjglse Detec ti on in thg R an dom Reg ion 
In an efficient TDMA system, the number of random 
(information) bits in a frame far exceeds the number of 
fixed bits. The majority of terras in the expression fer the 
probability of one or more false detections per frame. 
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eguation C.10, are .therefore of the form given by eguation 
C.4, the probability of false detection when only random 
bits are in the detector. With h=2. 



-3 

P{A(tH = 3.693 x 10 for N = 15 
r 

-4 

= 2.012 x 10 for N=20 



2. C ro ss- Co rrelat io n Bo upft 

When one or more unique uord bits, other than the 
one for which the detector is set, are contained in the 
detector, the cross-correlation criteria establish an upper 
bound on the probability of false detection (equations 3.4, 
3.6 and 3.7). These individual event probabilities can in 
turn be inserted into equation 3.8 to obtain an upper bound 
for the probability pf one or more false detections ( F * > as 
any one unigue word passes through the detector. Several 
such calculations were performed; for N=15 and N=20, the 
results were. 



-3 

F* = 5.831 x 10 for N=15 
max 

-3 

= 6.332 x 10 for N=20 



Thus F • is uithin an order of magnitude of the 
individual event probabilities in the random region. 
Because of the typically large number of random events, F * 
in fact makes an insignificant contribution to the total F 
over the period of one frame. Given the fact that a 
cross-correlation false detection pulse is potentially 
periodic, an "insignificant" F* is precisely what is 
desired. 

3 . Ex a S ystem Calc ulat io ns 

In crder to compute the probability of one or more 
false detections over a frame, it is necessary to have 
available the specific frame composition of the operating 
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system in question. A system was hypothesized for the 

purpose of obtaining a general feeling for the relative 
magnitude of the different factors involved. The number of 
stations was set at three in order to simplify the 
calculations and permit selection of a set of unique words 
from those available at N=15. The following set was 

selected at random: 



( 1 ) 0000011001011 11 

(2) 00010C100011 1 11 

(3) 0001010011011 11 



This set p.ossesses 


the 


following 


correlation 


properties: 









MI C COR RELATION 



in 



aiu 

ill 131 



ELATION 

2111 

121=111 121=111 121=121 



0 


15 


15 


15 


-14 


0 


0 


0 


1 


9 


9 


7 


-13 


0 


0 


0 


2 


7 


7 


7 


-12 


0 


n 


0 


3 


7 


7 


7 


-11 


0 


0 


1 


4 


5 


6 


6 


-10 


0 


i 


2 


5 


4 


4 


5 


-9 


2 


2 


2 


6 


5 


3 


6 


-8 


4 


3 


3 


7 


5 


5 


2 


-7 


5 


4 


2 


8 


4 


4 


3 


-6 


5 


5 


5 


9 


2 


2 


3 


-5 


7 


4 


5 


10 


1 


1 


2 


-4 


6 


5 


3 


1 1 


0 


1 


1 


-3 


5 


8 


6 


12 


0 


0 


0 


-2 


7 


8 


10 


13 


0 


0 


0 


-1 


1 1 


6 


8 


14 


0 


0 


0 


0 


1 1 


12 


10 










+ 1 


9 


10 


10 










+2 


9 


8 


8 










+ 3 


9 


8 


8 










+ 4 


3 


6 


5 










+5 


3 


3 


6 










+ 6 


5 


6 


4 










+ 7 


3 


3 


5 










4-8 


4 


4 


4 










+ 9 


2 


3 


3 










+ 10 


2 


2 


1 










+ 11 


1 


1 


1 










+ 12 


0 


0 


0 










+ 13 


0 


0 


0 










+ 14 


0 


0 


0 




The 


hypothetical 


system 


is 


patterned 


aft' 


tional 


system described in 


Ref 


. 3. The 


tra ns 


is 


60 


Mbps 


(4~phase) , 


the frame length is 750 us 



the 



guard band is 24 bj.ts in width and there are 60 fixed bits 
in the preamble for carrier and bit timing recovery. The 
latter are arbitrarily assumed to consist of a sequence of 
30 1*s for carrier recovery, followed by a sequence cf 30 
alternating 0's and I's for bit timing recovery. Bits are 
clocked in from the left, with the detectors set tc the 
unigue wcrds as listed above. 

It is first necessary to determine the number of 
random versus fixed detection events per frame. The product 
of the transmission hit rate and frame length gives the 
number of bits in one frame, n=45000. The beginning of each 
station’s burst is cpmposed of 75 fixed value bits (recovery 
plus unigue word) . When one or more of these bits is in the 
detector, the detection event is, by definition, ’'fixed.*' 
Therefore, for each station's burst, there are 75 fixed 
detection events as each bit, from the first to the last, is 
decked intc the detector; plus an additional 15 fixed 
events as the last bit, and those immediately preceding it, 
are clocked out of the detector. The total number of fixed 
detection events per frame is thus egual to 270 (3x90) . The 
number of random events is 44730 (45000^270). 

With this number of random events, the upper bound 
on the probability of one or more false detections pet frame 
in the random region alone is (equation 3. 10) , 

r3 

F < rP {A (t) } = 44730 (3.693 x 10 ) 
r 

F < 165 

which is net a particularly enlightening result for a 
probability upper bound. It is indicative, however, of the 
fact that a false detection is almost certain during 
acquisition. By way of comparison, the Ref. 11 expression 
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for the probability of no false detection in the random 
region gives, 

44730 -72 

1 - F = [ 1 - P{A(t) } ] =1.3x10 

r 

Therefore, it can be concluded that many false detections 
will occur during the acquisition phase, and that, for the 
system design hypothesized, little can be done to improve 
the situation, short of using some type of acquisition gate. 
This is generally true for the typical prototype TDM A 
systems which have been developed. No attempt was made to 
compute the expected acquisition time, a subject treated 
extensively in Hef. 11. 

If an acquisition gate is employed, the detector is 
active fpr a maximum of 2B=48 detection cycles pricr to the 
receipt of the associated unique word. All of these are 
fixed detection events. At each cycle, the autocorrelation 
values can be added to the number of matches occurring in 
the recovery bits to obtain the number of total fixed 
matches, k. Equation C.8e can then be used to obtain an 
upper bound cn Fw Dping so, for detector number (2) , 

-11 

F < 3.996 x 10 
10 

"F > 2.503 x 10 frames - 217 days 

For detector number (3) , 

-4 

F < 5.993 x 10 
T > 1669 frames = 1.25 sec 

The mean time between false detections (T) is more 
than acceptable for both, indicating a high probability of 
achieving acquisition in just one pass of the unique word. 
Acquisition cf the reference unique word, assumed to be 
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number (1) here, must still he effected without benefit of 
gating, and thus remains a potentially time consuming 



process. 



V. CONCLU SIO NS 



This research effort was undertaken for the express 
purpose of finding unique word sequences which would be 
suitable for use in a TDKA satellite communication system. 
The correlation properties of potential unique words were of 
primary concern in this effort, criteria being established 
which would insure a low probability of false detection. 
Numerous sets of unique words were found using these 
criteria . 

The results obtained can not, unfortunately, be applied 
to any given operational TDMA system without some additional 
consideration of false detection probabilities. In an 
attempt to generalize findings, it was necessary to ignore 
the fact that each s.ystem has its own peculiar set of fixed 
value bits, in addition to the unique words, which are 
transmitted in each station's burst. These fixed value bits 
are in fact an extension of the unique word sequence and 
produce a cross-correlation distribution of their okn when 
passing through a detector. The cross-correlation may very 
well be such that a set of unigue words selected at random 
from this work would prove unacceptable. Fortunately, the 
number of sets available of any given length, together with 
the existence of ;the complementary options, provides a 
flexibility sufficient to cover any foreseeable situation. 

The length of the sets found fell short of expectations. 
To begin with, in setting cross-correlation criteria based 
on worst-case possibilities for the number of detector 
agreements, unique words which might actually be useable in 
a particular systeu were excluded. That is, the number of 
agreements produced when non-unique word fixed bits are 
stepping through the detector can seldom be expected to 
approach the worst-case limit. This fact was noted at the 
outset, however, and accepted in order to generalize the 
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results as much as possible. The real disappointing factor 
was that, even with the strict correlation criteria, a 
computer search could not be designed tc uncover with 
certainty the sets of maximal length for values of N in 
excess of sixteen. 

There is a need for an efficient algorithm for combining 
unigue wcrds into sets possessing good cross-correlation 
properties. As stated in Kef. 11, the problem is a puzzling 
one-. It is a problem which can undoubtedly be solved only 
by adopting some type of short-cut technique, wherein but a 
few of all the multitude of different matching combinations 
are considered. The ability to predict the maximal length 
for a given N would be of tremendous benefit in this regard. 

The ultimate achievement would be the discovery of an 
underlying structure which enabled one to extract unigue 
word seguences without the need to resort to brute force 
computer search techniques. 
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APPENDIX A 
UNI2UE WORD SETS 



Selected sets pf unique words, representative cf the 
many found, are presented herein. The listing for N=15 is 
complete. In the interest of space, unique words are listed 
in octal format. All sets are valid for a detector 
threshold setting of two. 

A “gated set" is one in which the first sequence listed 
(the supposed reference) cross-correlates well with all the 
sequences which follow. A “complete set" is one in which 
each seguence crpss-correlates well with every other 
seguence . 



N = 1 5 



MAXIMAL LENGTH GATED GET 



01457 01347 02067 02327 02367 02467 03067 

04437 05157 04157 04657 04517 04117 



COMPLETE SETS 



(1) 00627 02067 02217 

(2) 06627 02217 02327 

(3) 06627 02217 02367 

(4) 01317 01667 02327 

(5) 01457 04437 05157 

(6) 01457 04437 04157 

(7) 01457 04437 041 17 

(8) 01667 02327 02637 



(9) 01667 02367 02717 

(10) 01667 02367 02637 

(11) 01667 03367 02637 

(12) 05477 04757 03557 

(13) 05477 03557 05157 

(14) 05477 03557 C4157 

(15) 03137 04437 05157 

(16) 04437 04157 03057 
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MAXIMAL LENGTH GATED SET 



002637 



001627 
002367 
003327 
004167 
004657 
006467 
014337 
01 1637 
010757 
010157 
010457 



001647 

002467 

003347 

004327 

004717 

006567 

010337 

011237 

012357 

011657 

010717 



N = 1 € 



001667 

002647 

003367 

004367 

004727 

006767 

011537 

014237 

010357 

006657 

011517 



002247 

003067 

003557 

004517 

005167 

014277 

014537 

010237 

006557 

010657 

010517 



002327 

003127 

003667 

004617 

006167 

011737 

006137 

006757 

005157 

006257 

010617 



002347 
003267 
004157 
004627 
006267 
012337 
014137 
004757 
006157 
01 1457 



COMPLETE SETS 



(1) 


001457 


0 0 2 J 1 7 


005317 


011237 


010657 


011517 


(2) 


001667 


002637 


004657 


005161 


006467 


006157 


(3) 


002367 


005.167 


006467 


005477 


006137 


006257 


(«) 


002467 


002637 


003267 


004367 


006137 


006257 


(5) 


002637 


003267 


004367 


006467 


006137 


006257 


(6) 


002637 


003267 


005167 


006467 


006137 


006257 


(7) 


C02637 


004657 


006137 


011237 


010357 


011457 


(8) 


002717 


005317 


010337 


003137 


010657 


011517 


(9) 


003057 


004567 


005067 


006267 


004437 


005157 


(10) 


003267 


004367 


004717 


005467 


003137 


005157 


(11) 


003557 


011537 


005637 


010557 


0 11657 


011317 


(12) 


004567 


005067 


006267 


006137 


004437 


005157 


(13) 


004417 


003137 


004437 


005157 


010157 


011517 


(14) 


005317 


011077 


010337 


003137 


010657 


011517 
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N = 1 7 



COMPLETE SETS 



(1) 


001667 

014457 


006257 

012157 


011057 

011437 


011167 
0 14237 


012467 


(2) 


003137 

011767 


005567 

012667 


005717 

012717 


007157 

013167 


007327 


(3) 


00332 7 
012667 


005567 

012717 


007157 

013167 


007267 

021017 


011767 


<4) 


002637 

024.157 


021137 

022257 


014137 

014257 


011437 

021457 


024437 


(3) 


002717 

007317 


003237 

011727 


005567 

012267 


006327 

013167 


007267 


. (6) 


.0031 37 
011767 


005667 

012267 


005717 

012717 


006727 

013167 


007157 


(7) 


005567 

011727 


005717 
0 12667 


006657 

012717 


007157 

013167 


007267 


(8) 


00571 7 
011767 


006657 

012667 


006727 

012717 


007267 

013167 


007567 









N=J6 ' 






COMPLETE SETS 


(1) 


001457 

020567 


004767 

021667 


01 1637 
025517 


0 14657 
014337 


015527 

024637 


(2) 


004657 

016667 


004767 

021767 


005717 

026317 


012717 

006277 


015727 

014337 


<3} 


012717 

006277 


014657 

014337 


025517 

045437 


044177 

041357 


040677 

022357 


<<») 


014657 

053637 


030277 

045437 


042477 

042757 


043337 

022357 


050337 

044357 


(5) 


021557 

050637 


044177 

045437 


042677 

041357 


024677 

022357 


030537 

031257 
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N = 19 



COMPLETE SETS 



(1) 


0004657 

0010767 


0005567 

0011167 


0005637 

0030577 


0006767 

0042277 


0007557 

0022477 


(2) 


0004657 

0104577 


0005637 

0306677 


0007557 

0014277 


0013357 

0101477 


0015757 

0022477 


(3) 


0007557 

0104577 


0011537 

0106677 


0013357 

0014277 


0013477 

0101477 


0015757 

0022477 


(4) 


0007557 

0104577 


0013357 

0101177 


0013477 

0106677 


0015757 

0014277 


0046577 

0101477 









N =20 






COMPLETE 


SETS 










(1) 


0004437 

0061177 


0006277 

0054677 


0023477 

0024677 


0030737 

0C43277 


0064577 

0040477 


(2) 


0004437 

0021137 


0007 137 
0024337 


0011537 

0C30737 


A ^ ~> o *3 *7 

v/ vy i o z. j # 

0030277 


0015637 

0045477 


(3) 


0007137 

0032737 


0011077 

0060577 


0013237 

0054677 


0015637 

0047277 


0021337 

0041477 


(4) 


0013237 

0035737 


0013477 

0-062577 


0014277 

0050677 


0020677 

0042277 


0032737 

0047477 


(5) 


0021137 

0064377 


0023077 

0035577 


0024477 

0054577 


0024737 

0036677 


0030337 

0043477 
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APPENDIX B 

Determination of bbrst d rift rate 



In the absence of transmitter burst position corrections 
and under zero doppler conditions, a given unique word will 
be received every T seconds. (Other factors, such as clock 
stability and logic jitter, affect this periodicity; 
however, they are inconsequential compared to potential 
doppler effects and will be ignored.) Burst drift is 
defined as the doppler induced, real-time shift of a burst 
irom its periodic position. 

Burst drift rate is a function of the relative velocity 
between transmitting station, satellite and receiving 
station. Interest here is in determination of the maximum 
drift rate that can be expected. First, let us consider the 
motion of the typical communication satellite in an 
equatorial, geosyncnjronous orbit. Ideally, such a satellite 



remains at the same altitude over a given point on the 
equator. In reality.-, it is physically impossible to achieve 
and maintain a perfectly synchronous orbit [12]. Reference 
13 derives equations which can be used to estimate the 
worst— case relative velocity between a geosynchronous 
satellite and a given point on the earth's surface. The 
orbj-t eccentricity and inclination are the two primary 
factors which influence the satellite's motion. 

A non-circular orbit results in a sinusoidal variation 
of the satellite altitude. To a first order approximation, 
the resulting doppler shift is independent of earth station 
coordinates. If e is the eccentricity of the orbit, the 
maximum relative vel.ocity is approximately. 



v = 12100e Km/hr = 3361e m/sec 
e 
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Using .0004 frpm Table 1 of Ref. 14 as a typical 
eccentricity for later generation communication satellites, 

v =1.3 m/sec 
e 

When the orbit is inclined from the eguatorial plane, 
the satellite moves in a figure-8 pattern, elongated north 
and south. The nprth-south motion predominates and is 
sinusoidal in nature. If i is the orbit inclination in 
radians, the resulting maximum relative velocity as seen by 
a station at latitude L can be approximated by, 

v = 2092isin(L) Km/hr = 581isin(L) m/sec 
i 

Satellite visibility extends to approximately 80° 
latitude. At this latitude, the satellite appears just 
above the horizon and relative velocity due to orbit 
inclination is a maximum. (It is in fact highly unlikely 
that a station would attempt to use a satellite below 10° in 
elevation.) Thus, assuming a maximum inclination cf 3° 
(.052 rad) and L=80°« 

v = 29.7 m/sec 
i 

A pessimistic estimate of the maximum relative velocity 
due to satellite motion alone can be obtained by combining 
the effects of orbit eccentricity and inclination, giving, 

.v = v + v = 31 m/sec 
s e i 

Next, the contribution of earth station motion must be 
considered. Reference 13 analyzes the headway, roll and 
pitch components of ship motion and concludes that the 
maximum possible relative velocity contribution is, 

v =38 m/sec 

g 
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For a shipboard platform, then, the maximum relative 
velocity between ground station and satellite is, 

v - v +v = 69 m/sec 
g s 

With representative numbers for satellite and surface 
platform motion in hand, it is next necessary to consider 
how they might combine to produce a shift in the position of 
a received unique word. The case of a station receiving its 
own unigue word (transmitted up to the satellite and back) 
is a simple tuo-statjon relative motion problem and the 
well-known dcppler eguation applies. 

d = 2vf/c (B.5) 



where f is the transmitted bit rate and c the speed of 
light. 

The case of a station receiving another's unique word is 
a three-station motion problem. Maximum doppler occurs when 
the maximum relative velocity exists between both ground 
stations and the satellite simultaneously; that is, when 
maximum up-link doppler is combined with maximum down-link 
doppler. Such a situation is produced when both ground 
stations are at the same 80° latitude and both are either 
opening or closing the satellite at the maximum rate. For 
this case also, then, with maximum relative velocity, v, 
existing between b oth stations and the satellite, eguation 
B.5 applies. 
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APPENDIX C 

DET ECTION PROBABILITIES 



Expressions for the probability of missed and false 

detections are derived in the following paragraphs. These 
expressions depend somewhat on general detector operation, 
but not to the extent that their application is restricted 
narrowly to any one particular system. 

A. HISS-DETECTION PROBABILITY 

If the bit error rate (BER) , i.e„ the probability that a 

received bit will be in error, is p and the number of bits 

in a unigue word is N, then, assuming errors occur 

independently, the probability of a unigue word being 

N 

received error-free ,is (1-p) . With the detector threshold 

set at h, a unigue word will be detected if it contains h or 

fewer errors. The probability of detection, P, is thus 
obtained by summing terms in the binomial distribution which 
correspond tc the probability cf no errors, the probability 
of one error, up to the probability of h errors. Using 

C(N,i) to represent the binomial coefficient N!/i!(N-i)i, P 
is given by, 

h i N-i 

P = ? C(N,i)p (1-p) (C.1) 

1=0 

The miss-detection probability (the probability of 
receiving a unigue word with more than h errors) is given 
^Y/ 

N i N-i 

Q = 1-P = 2 C (N,i) p (1-p) (C . 2) 

l=h+ 1 

The mean time between miss-detections is equal to the 
reciprocal of the expected value of "number of misses per 
unit time." As there is only one possible miss-detection 
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event per frame, .the expected value is equal tc the 
“iss- detection probability, Q, and the mean time between 
miss-detecticns is. 



-n _1 “I 

M = Q frames = TQ secs 



(C.3) 



B. PI20BAEILITY OF FALSE DETECTION 

The following .notation will be used in developing 

expressions for false detection probabilities: 

'A (t) 'A (false detection event at time t. 

Used interchangeably with "A . " 

t 

F {A ( t ) } The probability of A (t) . 

E (t) No false detection at time t. 

Used interchangeably with B 

t 

P{B(t)} The probability of B (t) . 

t An integral time variable, the basic unit 

of which is one bit cycle (t=ib) . 



On any given bit cycle, other than the one in which the 
actual unigue word is completely contained in the detector, 
a false detection pulse will be produced if N-h or mere of 
the N bits in an active detector match with the detector 
setting (B>h) . The probability of this occurring depends on 
whether or net the value of transmitted bits are fixed or 
variable. The beginning portion of each burst is composed 
of fixed-value bits (timing recovery and unigue word) , while 
the remainder of the burst contains bits of variable value. 
The probability of one of the latter being transmitted as a 
•0* or a • 1 • depends on the type of information transmitted. 
As a first order approximation, it will be assumed that the 
variable bits are in fact random, taking on the value *0* 
and 1 1 * with egual probability. With random bits alone in 
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the detector, the probability of a false detection cn any 
given bit cycle is, 

P{A(tn =.£** C(N,i) (1/2) 1 ( 1/2) N 
r i=0 

N h 

= (1/2) 5 C(N,i) (C.4) 

1=0 

The probability of false detection is not so easily 
computed when there are fixed value bits in the detector. 
"Fixed value" is used to denote those bits in the preamble 
which are transmitted with the same value {*0' or 1 1') in 
each frame. Because of the presence of transmission errors, 
these bits may actually be received with the opposite value, 
the probability of which is p. An expression fcr the 
probability of false detection can be derived based on the 
number of fixed bits in the detector and the number that 
matched with the detector setting when transmitted. Such an 
expression must account for all the possible combinations of 
errors within the fixed bits and as such will include 
additive terms containing higher powers of p. With p<<1, a 
good approximation can be obtained by discarding all terms 
which contain other than the lowest power of p. Therefore, 
rather than derive a complex general expression from which 
an approximation can be drawn, the condition p<<1 will be 
invoked at the outset in order to simplify the derivation. 

Suppose there are m fixed bits in the detector and that 
k of them matched with the detector setting when 
transmitted. A false detection will occur when the number 
of agreements within the fixed bits, plus the number of 
agreements within the remaining N—m random bits is greater 
than or egual to N-h,. With p<<1, interest can be focused on 
those particular combinations which require the fewest 
number of transmission errors in order to produce a false 
detection. Three separate cases will be considered. 
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CASE 1. The number of fixed agreements, k, is 
greater than or egual to N-h. A false detection is assured 
if the m fixed bits are received error free, and the 
probability cf false detection can be approximated by, 

P-{A (t) J ^ = ( 1 — p) k + h > N (C.6a) 

CASE 2. The number of fixed agreements, k, is less 
than N-h, but enough random bits (N-m of them) remain to 
produce a total of N-h or more agreements. That is, 

k+N-m>N-h 



A false detection will be produced when there are no errors 
in the fixed bits and N-k-h or more agreements within the 
random bits; or, equivalently , when there are 

N-m - (N-k-h) = k + h- m 

or fewer disagreements within the random bits. The 
probability of false detection can thus be approximated by 
multiplying the probability of receiving the m fixed bits 
error free by the probability that the N-m random bits will 
contain k+h-m or fewer disagreements. 



P{A (t) } = 

£ 

w N-m k+h-m 

(1-p) (1/2) 2 C (N-m , i) ; m<k+h<N (C.Eb) 

i=0 



CASE 3. The number of fixed agreements, k, plus the 
number cf random bits, N-m, is less than N-h; i.e., k+h<m. 
A false detection can not occur unless there are errors in 
the fixed tits. The fewest number of errors are reguired 
when all cf the random bits and all of the original k bits 
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match with the detector setting. In this case, enough 
errors must cccur in the remaining m-k fixed bits to produce 
N-h agreements (the minimum required for false detection) . 
The number of errors required is, 

a = N - h - (k + N-m) 

= m - k - h 

The probability of false detection can thus be approximated 
by taking the product of (1) the probability that the N-m 
random bits will all match with the detector setting, (2) 
the probability of receiving the original k matching bits 
error free and (3) the probability that the m-k fixed bits 
will contain 'a' errors. 



P{A(t)} £ 



U-m k a ro-k-'a 

= (1/2) (1~p) C (m-k, a) p (1-p) 



N-m a ra-a 

(1/2) C (m-k , a) p (1-p) 

N-m m-k-h k+h 

(1/2) C (m-k, ra-k-h) p (1“p) 



k + h<m 
(C. 6c) 



If all the bits in the detector are fixed, the 
probability of false detection in any given bit cycle is 
given by eguatior^s C.8a and C.8c with m=N. 



N 

P { A (t) } = (l-p) ; k+h>N 



(c. ed) 



P(A(t) j 



N-k-h k+h 

c (N-k, N-k-h) P (1-p) ; 



k+h<N (C„8e) 



1 * Joint Probabil ities 

Equations CU 4 and C.8 give the probability of false 
detection for a single detection event within the frame. 
Calculation of the probability of false detection over 
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several consecutive bit periods, specifically over one 
frame, is complicated by the fact that successive detection 
events are net independent; i.e. , 

P{A(t) ,A(t+1))*P[A(t)}P[A (t + 1)} . 

That this is true can be seen by noting that successive 

detection events are generated by shifting the contents of 

the shift register by one bit. If a false detection does 

(or does net) occur at one step, possibilities for the 

contents of the shift register are necessarily restricted to 

certain combinations of bit values. These bit values carry 

over to the next detection event and reduce the sample space 

accordingly. As a sample, extreme example, consider the 

case of the urjigue word 101010. The probability of false 

detection with random bits in the detector and a zero 

6 

threshold setting is, from equation C.4, (0.5) . However, 
given that a false detection occurs (the detector contains 
101010), the detector will contain X10101 on the next step 
and the probability of false detection is zero. 

Dependent joint probabilities are best dealt with by 
locking for, and finding, grounds for the assumption of 
independence. Schrempp and Sekimoto [11] based their 
independence assumption on the previous uork of Hill [14]. 
The validity of this assumption is nearly impossible to 
access without making lengthy comparative calculations for 
unigue words of different length, a task of monumental 
proportions. Bather than attempt to justify such an 
assumption here, jpint probabilities will be expressed in 
terms of an upper bound. 

There are n detection events per frame during 
acquisition, n-1 of which may produce false detections. 
Defining t=0 at the point when the associated unique word is 
completely contained in the detector, the probability of no 
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false detection during a frame period is given by. 



H = P{E -B ,***,B } ( 3 - 9 ) 

12 n-1 

= P{E }P{B /B } P [B /B B }*«<>P{B /B B •••B ) 

1 2 1 3 12 n-1 .1 2 n -2 

Cr, from a different viewpoint, if F is the 

probability cf one cr more false detections per frame, 

H = 1 - F 



n-1 n-1 n-1 

* -. 2 - P(A.) “ .2, 2. P(A. ,A .) 

i=1 l 1=1 3=1 l 3 

n-1 n-1 n-1 

+ . 2 . . 2 . Jl . P (A , , A , A ) - 

i=1 3=1 k=3 13k 

± P (A , A , • ® • , A ) 

12 n-1 



For straight-forward calculations, F is obviously 
orders of magnitude more difficult to obtain than is ft. 
However, its form does offer possibilities for reducing the 
complexity cf .the problem. Although successive events are 
not independent, events separated by more than N detection 
cycles are - - after N detection cycles, all bits from the 
first event have been shifted out of the detector. For 
j>i+N , 

P {A , A ) = P {A } P {A } 
i 3 i 3 

and similarly for higher order joint probabilities. If 
P{A(i)}<< 1 , these product terms are small compared to the 
individual terms from the first summation in F. However, 
for the dependent .joint probabilities it is guite possible 
that 



P { A , A } = P{A }P{A /A } 

i 3 i j i 
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is of the same order as P {A } . The same can be said for 

i 

higher crder terms, but with a decreasing expectation cf its 

actually proving true. It is therefore not possible to 
discard as insignificant all of the terms involving joint 
probabilities. The most that can be said is that, 

P{A i>P{A , A }>P{A , A , A } . . • 

1 i j i j k 

And that. 



F 



n- 1 




P {A . } 
x 



The false detection probabilities P { A (i ) } are 
obtainable from equations C.4 and C.8 as appropriate. For 
the random detection events (equation C.4) , 

P {A (i) } =t {A ( j) ] ; for the fixed detection events, the 
probabilities must be computed individually. If there are r 
random false detection events and f~n- !-r fixed false 
detection events per frame. 



F 

ma x 



m : 1 



P{A.} 

x 



= rP{A(t)} + £P{A(t)l (C. 10) 
r r f 

When the detector is gated during acquisition, as 
described in Section II. C, the detector is active for only 
2B+1 bit cycles about the nominal position of the associated 
unigue nerd. The true detection pulse may occur anywhere 
within the aperture, depending upon the amount of 
displacemen t of the unique word from its nominal position, 
and the maximum number of false detection events is equal to 
2B. The probability of one or more false detections in the 
aperture prior tc the receipt of the unigue word is 
therefore of the same form as equation C.10, with the number 
of events equal to 2B vice n-1. 
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In the steady-state mode, the aperture width is 2D+1 
and the maximum number of false detection events is equal to 
2D. 

In summary, equation C.10 provides a bound on the 
probability of one or more false detections over one frame 
period, or, with a change of limits, over an acquisition or 
steady-state detection aperture. This bound is useful in 
estimating the frequency of false detections and in 
assessing the relative merits of different unique word 
combinations . 

As an aside, it is interesting to compare the random 
event term of equation C.10 with the results of Ref. 1-4, 
obtained by assuming independent events. For random bit 
detection events, it was concluded that, 

F = 1 - [1 - P £A (t) } ] r 



Expanding, 



2 

F = 1 - [1 - rP {A (t) } + C (r , 2) P {A (t) } - •••} 



For F { A (t) } « 1 , a basic assumption of Ref. "\ii, an 
approximate expression for F is, 

F = rP {A (t) } 

Which is identical to the random bit term in equation C.10. 

2. Mean Tim e get ween. Fal se De tections 

A false detection, both in the acquisition and 
steady^state modes, deactivates the detector and sets an 
aperture for the next frame. The probability of cycling 
through several successive detection events without a ny 
false detections is therefore of primary concern. Fron the 
opposite viewpoint, it is the probability of one or more 
false detections over a given time span which is meaningful, 
not the probability of false detection at any given step. 
This distinction is important in deciding upon what tc use 



as an appropriate measure of "mean time between false 
detections" (MTBFD) , and in turn interpreting its meaning. 

If the detector were not deactivated by a false 
detection pulse, several false detections might occur over 
the period of one frame. The MTBFD would then be equal to 
the reciprocal of the expected value for "number of false 
detections per framei" With detector deactivation, only one 
false detection per frame can actually occur and the MTBFD 
is affected accordingly. In this case, for computational 
purposes, the probability of a real, physical false 
detection pulse occurring somewhere within the frame is 
equal to F, the probability of one or more false detections 
per frame. The MTBFD is therefore equal to the reciprocal 
of F and can be bounded by substitution from equation C.10. 



a? 



F > F frames 
max 



(C. 11} 
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APPENDIX £ 

COMPU TER f ROGRAf} 

The computer program was written, and debugged, in three 
separate functional segments and then patched together to 
fgrm one continuous program. These segments will be 
described and listed separately. The last segment listed is 
the basic cross-correlation matching routine. Several 
different revisions were employed in an attempt to overcome 
time and memory limitations. 



A rra ys 



GUW (IEIM) - Sequences which pass the autocorrelation 
test . 

GCTUW (IDIM) - Complement transpose of GUW. 

GTUW(IDIK) - Transpose of GUW. 

GCUW (IDIM) - Complement of GUW. 

MTEME(IDIM2) - Used for temporary storage of good 
crcss-correlation matches. 

INCR ( ) - Used to store pointers needed in the 

crcss-correlation matching process. Dimension greater 
than or equal io maximal length set expected. 

MATCH (IDIMl,IDIM2) - Matrix in which good 
crcss-correlation matches are recorded. 

MCROSS (2IDIM) - Sequences from autocorrelation sort 
which are checked for cross-correlation (the GUW and 
GCTUW subsets},. 



Equivalence Statement - The first half of MCROSS is 
stored in same memory location as GUW( ); the second 
half in same location as GCTUW ( ) . 
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V ar i ah les Co a a on to All Segm ents 

UW - Seguence under test. 

SUW - Shifted version of UW. A seguence is shifted one 
bit at a time and, at each step, added mod-2 to UW. 

TUW - Transpose of UW. 

CTUW - Complement transpose of UW. 

CUW - Complement of UW. 

N — Length of sequences being tested (Input) . 

IEPS - Threshold.; h of thesis main body (Input) . 

IALFA ~ Aperture parameter; equals 2D (Input) . 

NUW - Number of good UW's obtained in autocorrelation 
sort. 

NAUTO - Number of UK's and CTUW's, combined, which are 
checked for crpss-correlation; equals 2NUW. 



67 



. Auto co rr elation Sc^rt 



Var iab le s 



K - Smallest integer equal to or greater than N/2. 

LAST,LAS11 - Intermediate skip points in the sequential 
checking process. Sequences which must be checked 
tend to occur in blocked subsets, with two redundant 
seguences separating each. 

IZ1 - Next to last sequence which must be checked. 

IZ2 - Last sequence which must be checked to insure all 
possible non-redundant sequences have been considered. 

MASK - Binary constant consisting of N I's, right 

justified. Used to obtain the complement of an N bit 
seguence. 

IFLAG - Set to one when last sequence being checked. 

Keys the termination of autocorrelation search. 

JFLAG - Set to one when autocorrelation condition 

reguires R(t)=Q. Autocorrelation computation for a 
particular UW is complete when JFLAG=IEPS. 

I - Shift variable. Equivalent to time variable of 
thesis main body. 

ICOUNI - The autocorrelation value at step I. 



Ou tp ut 



The seguences which pass 
the corresponding CUN's, TUW's 
side-by-side decimal and octal 



the autocorrelation test, 
and CTUW's, are output 
format. 



and 

in 
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CO <0 (A CO CO CO CO (rt CO CO CO CO u 






EGER UW,SUW,TUW,CTUW,CUW,GUW (73) # GCTUW(73) , 
gj73)^GCUW (. 73 ) , MTEMP (110) ,INCR ( 36 ) , MATCH ( 76 , 1 10) , 

(1) , MCROSS (1) ) , 

PARAMATER iNPU'i 



, (GCTUH (1) 



MCROSS (74) ) 



OF NEXT UW 



INTEGER 
1GTU 
2MCEO 

EQUIVALENCE (GUW (1) , MCROSS (1 

C INITIALIZATION AND A 

IDIM = 7 3 
IDIM 1=73 
IDIM2= 110 

1 READ (5,500) N, IEPS , IALFA 

WRITE (6,505) N, IEPS, IALFA 
IF (N.EQ. (N/2) *2) GO TO 5 
•K— ( N +' 1 } / 2 
LAST=2**K-2 
LA ST 1 = 2** (K+1) -2 
GO TO 10 
5 K= N/2 

LAS1=2** (K+1) -2 
LA ST 1 = 2** (K + 2) -2 
10 IZ1=1 

MLAST= N-2 
DO 15 M=K, MLAST 
15 IZ 1=1 Z 1+2**M 

IZ2=2** (N-1) -2 
MASK=2**N— 1 
IFI AG= 0 
UW = 0 
3= 1 
L= 1 

J r L AG- 0 

C INITIALIZE FOR AUTC-C 
20 1=0 

SU K=U H 

C COMMENCE AOTO-C LOOP. 

C THE FOLLOWING IN-LINE SEGMENT INITIALIZES AN N-I COUNTER, 
C LOADS SUW INTO A -REGISTER , SHIFTS SUN RIGHT ONE BIT, 

C STORES RESULTS IN "SUW" FOR LATER USE, COMPUTES MOD-2 SUM 
C OF UW AND SUW, COUNTS NUMBER OF ZEROS IN FIRST N-I 
C POSITIONS (FROM RIGHT) OE THE MOD-2 SUM. 

30 1=1+1 
ICCUNT=0 
LD A N 
SUE I 

ADD =077677777 
RCH 070100,1 
IDA SUW 
SHIFT C21001 
STA SUW 
EOfi UW 
SHIFT 0100-1 
SKE =040000000 
MPC ICCUNT 
BEX $-3,1 

ARE WE INSIDE APERTURE? 

IF (I. LE. IALFA) GO TO 31 
C ARE WE ON DISTRIBUTION TAIL? 

IF (I.GE.N-IEPS) GO TO 32 

C DOES UK PASS AUTO-C TEST AT STEP I OUTSIDE APERTURE? 

IF (ICCUNT. LT.N-IEPS-I) GO TO 30 
GO TO 35 

C DOES UK PASS AUTO-C TEST AT STEP I INSIDE APERTURE? 

31 IF (ICCUNT. LT.N-2*IEPS-I) GO TO 30 
GO TO 35 

C DOES UW FAIL ZERO TEST ON DISTRIBUTION TAIL? 

32 If (ICCUNT. NE.O) GO TO 35 
C NOW ON DISTRIBUTION TAIL 
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33 JFLAG= JFJ.AG+ 1 

IF (JPLAG.LT. IEPS) GO TO 30 
C HAVE G CCD UW. COMPUTE CUW, TUW, CTUW. 

C THE FOLLOWING IN-LINE SEGMENT INITIALIZES AN N COUNTER, 

C LOADS UW INTO A-REGISTER, FORMS TRANSPOSE OF UW BY RIGHT 
C SHIFTING ONE BIT AT A TIME INTO B-REGISTER AND THEN CYCLE 
C (END-ARCUNL) SHIFTING THE LATTER. FORMS COMPLEMENT OF 
C BOTH UW AND TUW USING MASK OF ONES. 

S LD A N 

S ADD =077677777 

S RCH 070100,1 

S LD A UW 

S RCH 037710 

S SHIFT C1001 

S SHIFT C 16002 

S BRX $-2,1 

S SHIFT Cl 1001 

S STE TUW 

S LDA UW 

S EOR MASK 

S SI A CUW 

S LDA TUW 

S EOR MASK 

S SI A CTUW 

C HAVE WE PREVIOUSLY CHECKED AND STORED THESE SEQUENCES? 

IF (TUW. LI.UW.Ofi.CTUW.LT. Uri) GO TO 35 
C DIMENSION OVERFLOW CHECK 

IF (L.G3C.IDIH) GO TO 4 7 
C STORE GCOD SEQUENCES 
GUN (L) =U w 
GCUW (L) =CUW 
GTUW (L) =TU W 
GCIUW (L) =CTU W 
L=L + 1 



C PICK A NEW SEQUENCE 
35 IF (Uw .LT.LA5T) GO TO 40 

IF (UW.GE.IZ1) GO TO 45 
IF (UW . IT . LAST 1) GO TO 52 
J = J + 1 

IF (N. EQ. (N/2) *2) GO TO 50 
LAST1=2** (K+J) -2 
GO TO 52 

50 LAST1 = 2** (K+J+1)-2 

52 LASl=LAST+2** (K-J) 

UW=UW+3 

55 GO TO 20 

40 UW=UW+1 

GO TO 20 

45 IF (IFLAG.EQ. 1) GO TO 60 

IFL AG=IFLAG+ 1 
U W = IZ2 
GO TO 20 

C DIMENSION OVERFLOW HALT 
47 WRITE (6,506) UW 

STOP 

C COMPUTE NUW AND NAUTO 
60 L 1=1- 1 

N U W = L 1 



N A UTO= 2*N UW 

C OUTPUT GOOD SEQUENCES (DECIMAL AND OCTAL FORMAT) 

WRITE (6,507) 

65 WRITE (67568) M,GUW(M) ,GUW(M) ,GCUW(M) ,GCUW(M) , GTU K (H) , 

1GTUW (M) , GCTU W (M) ,GCTUrf (M) 



C 

500 

505 

506 

507 

508 



FORMAT N = ' ,12,1 OX, 'IEPS=' , 1 1 , 1 0 X , ' IALFA= ' ,12) 

FORMAT ( ' 0* » ' AUTOCORRELATION SEARCH TERMINATED TO , 
i' PREVENT 6lMENSION OVERFLOW. ; LAST WORD CHECKED. ,-6) 
F O I? N A T MU* 1 UNIQUE WORDS W H J. C U PASS AUTOCORRELATION r 
I 'TEST V/l6x< U^' /25X, ! CUW 1 , 2 5 X , 'TUW , 2 3 X , 1 CTUW'//) 
FORMAT {' ' ,l6,5X,4 (16, 3X, 010,8 X) ) 
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Cross- Correl at ion Sort 



Va riable s 



M2 - Number of sequences with which sequence being 
tested has cross-correlated acceptably. Used for 
filling MATCH .matrix. 

M2HIGH - Highest number of sequences with which any one 
sequence has cross-correlated acceptably. Used for 
MATCH matrix output. (This is not the maximum number 
of cross-correlation matches existing. Correlation is 
checked, and matches scored but once, commencing with 
the first sequence in the AUTO-C set and working up. 
Thus M2HIGH does not include matches betueen the 
seguence being tested and previous sequences in the 
AU'IO-C set,) 

Other variables are defined under "Arrays" and 
"Comiscn. " 



Output 



(1) The UW and GTUW sequences from the AUTO-C set which 
are being cross-correlated are output in decimal format, 
together with an index number which denotes the order in 
which the sequences are tested. 

(2) The matching matrix (MATCH) is output in 
two-dimensional form. The elements of MATCH are the 
index numbers corresponding to the AUTO-C set of 
sequences. There is one row for each AUTO-C UK (See 
Section IV.B). The first element cf a row is the index 
number of the UH, and the remainder of the row contains 
(higher) index numbers of the U2*s and CTUW's with which 
it matches. 
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C)ininU)U)lninininu)tninin -*>non Citsiinlnininunnuiin 



P rog raa 



C FILL MATCH KITH ZEROS 
DO 103 I=1,IDIM1 
DO 103 J=1,IDIM2 

103 MAICH(I,J)=0 
C PACK MCROSS 

DO 104 1=1. NUW 

C OUTPUT THE AUTO-C SET ELEMENTS AND CORRESPONDING INDICES 

104 MCROSS (NUN+I) =GCTUW (I) 

WRITE (6,610) 

WRITE (6,611) (I, MCROSS (I) ,1=1 ,NAUTO) 

C INITIALIZATION 
L= 1 



M2 = 2 

M2HIGH=2 

C SELECT NEXT UK FROM AUTO-C SET 
105 UK = MCRCSS(L) 

C DIMENSION OVERFLOW TEST 

IF (L.GT.IDIM1) GO TO 155 
C STORE INDEX NUMBER OF UK 
HTEMP < 1) =L 

C SET INDEX OF FIRST AUTO-C ELEMENT BEING CHECKED AGAINST UK 
M=L + 1 



C CHECK CROSS- CORK EL ATI ON WITH ZERO SHIFT BETWEEN UN AND 
C SOW. THE IN-LINE SEGMENT IS IDENTICAL IN PURPOSE TO THE 
C FIRST IN-LINE ROUTINE OF PREVIOUS PROGRAM SEGMENT 
110 ICC UNT-0 
1=0 

JF1 AG=0 



SU i=MCROSS (M) 

LD A N 

tiU V -U / / Q/ / / / / 

RCH 070100,1 
LD A UK 
EOE SUW 
SHIFT 01001 
SKE =040000000 
MPC ICC U NT 
BRX $-3,1 

DCES PAIR FAIL TEST AT 1=0? 

IF (ICOUNT. GE.N-IEPS) GO 
COMMENCE LOOP FCR COMPUTING 
CORRELATION. THE IN-LINE S 
TC THE CUE ABOVE. 

15 1=1+1 



TO 1 17 

RIGHT-SHIFTED CROSS- 
EGMENT IS IDENTICAL IN PURPOSE 



ARE 



ICOUNI-O 
LD A N 
SUE I 

ADD =077677777 
RCH 070100, 1 
LD A SUN 
SHIFT C21001 
STA SUW 
EOE aw 
SHIFT 01001 
SKE =040000000 
MPC ICCUNT 
BRX $-3,1 
WE CN DISTRIBUTION 



IF (I. GE.N-IEPS) GO 
C DCES PAIR PASS TEST AT 

IF (ICOUNT. LT.N-IEPS-I) 
C increment INDEX 
117 M= M + 1 



TAIL? 

TO 140 
STEP I? 

GO TO 



115 
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C HAVE HE CHECKED UW AGAINST ALL HIGHER ORDER SEQUENCES? 

IF (M. LE„ NAUTO) GO TO 110 
C FINISHED HITH UW. STORE HATCHES. 

120 M2 = M2- 1 



DO 125 K=1,M2 
MATCH (I, Kj^ = MIEMP (K) 

_ 1< OF UH MATCHES IN FIRST 

MATCH (I,IDIM2) =M2 
M2 EIGHER THAN FOR ANY PREVIOUS 
IF (M2 . GT . M2 HIGH) M2 HIGH=M2 
M2 = 2 



125 MTEMPjK 
C STORE NUMBE 



C IS 



UNUSED ROW ELEMENT 
UW? 



C TIME TC GET NEW UW,. 

L=I + 1 

C ABE WE FINISHED? 

IF (L.LE.NUK) GO TO 105 
C CROSS-C SORT COMPLETE. OUTPUT MATCH. 
135 WRITE (6,655) 

DC^ 137 1=1, L 

137 WRITE (6,66 0) (MITCH (I , J) , J= 1 , M2HIGH) 
GO TO 300 



GC TO MATCHING ROUTINE 



40 



DCES PAIR PASS ZERO 



TEST AT STEP I? 
GO TO 117 



IF (ICO UNT. GT . Oi 
JFI AG=JFLAG+ 1 
ARE WE FINISHED WITH RIGHT-SHIFT? 

IF^UFLAG.LT.IEPS) GO TO 115 

JFI AG= 0 
SUW=MCSOSS (M) 

START LOOP FOR COMPUTING LEFT-SHIFTED CROSS-CORRELATION. 
THE IN-LINE SEGMENT INITIALIZES N-I AND I COUNTESS, IOADS 
SUW INIC A-RjiGISi’E-fi, SHIFTS 5UW LEFT ONE BIT, STORES THE 
RESULT ECR LATER UEE, COMPUTES MOD-2 SUM OF UW AND SUW, 
SHIFTS THE MOD-2 SUM EIGHT I BITS SO THAT THE N BITS OF 



INTEREST ARE RIGHT- JUSTIFIED 
NUMBER OF ZEROS IN FIRST N-I 
THE MOE-2 SUM. 



IN THE A-REGISTER, COUNTS 
POSITIONS (FROM RIGHT) CF 



45 



1=1 + 1 
ICC,UNT=0 
L DA N 



SUE 
ADD 
RCH 
LD A 
ADE 
RCH 
L DA 



=077677777 
070100, 1 
I 

=077677777 
070100,2 
SUW 

SHIFT 02 500 1 
STA SUW 
ECR UW 
SHIFT 01001 
BEX $-1,2 
SHIfT 01001 
SKE =040000000 
MFC ICC UNT 
BEX $-3,1 

IF (I.GE.N-XEPS) GO TO 150 
IF ICOUNT.LT. N-IEPS-I) GO TO 145 
GO TO 117 

150 IF (ICO UNT. GT . 0) GO TO 117 

JFI AG = JFLAG+ 1 

IF (JFIAG.LT. IEPS) GO TO 145 
C HAVE A GOOD MATCH. STORE IN MTEMP 
MTEMP (M2)=M 
M 2 = M 2 + 1 

C DIMENSION OVERFLOW CHECK 

IF (M2.IT.IDIH2-1) GO TO 117 
WRITE (6,665) L 
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cnO 



155 



GO TO 120 
WRITE (6,670) 
GO TO 135 



10 FORMAT (' 1' ,' AUTOCORRELATION WORDS FROM ABOVE WHICH', 

1' ARE EEING CHECKED FOR CROSS-CORRELATION MATCHES'//) 
611 FORMAT (6 (1X,I4,'. ' 2X,I5,8X) ) 

655 FORMAT j'O* 'GOOD CROSS-CORRELATION PAIRS. THE NUMBER', 
1' LISTED IS THE SUBSCRIPT OF THE CORRESPONDING GCOD', 
2' AUTOCORRELATION WORD.'//) 

660 FORMAT (' O' , 10 (IX, 2016/) ) 

665 FORMAT (' O' ,' OVERFLOW OF DIMENSION IDIM2 ATTEMPTED', 

1 ' WHEN L=' 3 1) 

670 FORMAT ('O' 'PROGRAM TERMINATED TO PREVENT OVERFLOW ' , 

1' CE DIMENSION IDIM1') 
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Match ing Routin e 



Variable s 



HOST - length of largest set found at any given point in 
search. Smaller sets, when found, are not output. 

LTES1 - Sets the last row element which must be looked 
at based on current value of MOST. 

MAT1 - An element from the matrix row which has M AT 2 as 
the leading element. 

MAT2 - The last element of the current set being formed. 

MAT3 - An element from the matrix row which has one of 
the set members as a leading element. Once it is 
known that a potential addition matches with the first 
and last elements of the set, MAT3 is used to check if 
it matches with the other elements as well. 

Other variables are defined under "Arrays'* and "Common." 



Output 



Match sets are output one line at a time in 
format with the autocorrelation index preceding each 
of the set. 



decimal 

element 
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C SETS ABE FORMED BY FINDING ALL THOSE THAT HAVE THE I«TH 
C BOW LEADING ELEMENT AS THE LEADING ELEMENT OF THE SET. 

C IT IS KNOWN THAT ALL ELEMENTS OF A ROW MATCH WITH THE 
C FIRST AND IT IS NECESSARY TO DETERMINE WHICH OF THESE 
C IN TURN MATCH WITH ONE ANOTHER. IF A ROW ELEMENT IS 
C GREATER THAN NUW, IT COMES FROM THE CTUW PORTION OF THE 
C AUTO-C SET AND WILL NOT APPEAR AS A BOW LEADING ELEMENT. 
C FOR THIS CASE, THE CORRESPONDING UW BOW IS USED TO 
C DETERMINE WHAT FURTHER MATCHES EXIST. 

C 

c 

C INITIALIZE 
300 WRITE (6 , 700) 

M0ST=3 



1= 1 
M ~ 2 

C START LOOKING FOR SETS WHICH HAVE AUTO-C NO. 1 AS LEADING 

C ELEMENT 

305 K= 2 

C COMPUTE I TEST 

310 IF (MATCH (I,K) .GT.NUW) , LTES T=MATCH (MATCH (I , K) -NOW, 
1IEIK2) +M+1-MOST; GO TO 312 
LTEST=MATCH (MATCH (I,K) ,IDIM2) + M+ 1-MOST 

C GO TO NEXT BOW 

312 IF (MATCH (I,K) .GT.O) GO TO 315 
1 = 1*1 

C ARE WE FINISHED? 

IF (I. GT.NUW) STOP 
GO TO 310 

C IS NEXT ELEMENT IN BOW A ZERO? IF SO, TIME TO GO TO 

C NEXT BOW. 

315 IF (MATCH (I, K+1) .GT.O) GO TO 320 
1=1 + 1 



IF (I. GT.NUW) STOP 
GO TO 305 
320 L= 1 

J=K+ 1 

C IS MATRIX ELEMENT A 'CTUN»? 

325 IE (MATCH (I, K) .GT.NUW) GOTO 

... II EXT 



C START LOCKING FO 



327 



L=L+ 1 
IF (L . GT . LTES T) 



ADDITION 



415 

TO MATCH 



SET 



GO TO 372 



330 



C DO 
335 



MAT 1 = MATCH (MATCH (I,K) ,L) 
* ~~ m '' 3 7 



TO 



370 

TO 335 



IF SO, TAKE THE 



IF (MAT 1. EQ„ 0) GO 
MAT2=MAICH (I ,J) 

IF (MAT2„LE„ 0) GO TO 
IF (MAT2.LE. MAT.1) GO 
GO TO 327 

WE HAVE A POTENTIAL ADDITION TO SET? 

IF (MAT 1 . EQo MAT.2) GO TO 340 
J=J + 1 
GO TO 330 

C IS THE CURRENT SET LENGTH EQUAL TO TWO? 

C NEW ADDITION WITHOUT QUESTION. 

C^MUST^CHECK^liPpOTENTIAL^NEW ADDITION MATCHES WITH OTHER 
C ELEMENTS ALREADY Hi SET 
M 1 = M- 1 

DO 350 MM=2 , Ml „ ^ - r _ 

IF (MTEKP (MM) .GT. NUW) GO TO 352 

'LL — 2 

MAT 3= MATCH (MTEMP (MM) , LL) 

FOLLOWING TRUE, N^W a5dITI0N NO GOOD 
IF (MAT3.EQ.. 0.O3.MAT3.GT.MAT2) GO TO 325 
IF ( M A T 3 - EQo M AT2 ) GO TO 350 
L L = I. L + 1 
GO TO 345 
350 CONTINUE 



345 
C IF 
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C HAVE NEW ADDITION TO SET 
GO TO 3 55 

352 M AT2=MAT2-NU W 

DO 354 MM1=MM,M1 
LL = 2 

353 MAT3=MATCH (MTEMP (MM 1 ) - NU W , LL) 

If (MAT3.EQ. O.OR.MAT3. GT.MAT2) GO TO 325 
IF (MAT3. £y„ MAT2) GO TO 354 
LL-LL+ 1 
GO TO 353 

354 CONTINUE 
MAT2=MAT2+NUW 

C STORE NEW ADDITION AND SET POINTER 

355 M= M + 1 

MTEMP (M) =M AT 2 
I NCR (M ) = J 

C START LOOKING FOR SEQUENCES WHICH MATCH WITH FIRST AND 
C (NEK) LAST ELEMENT OF SET. 

IF (MATCH (I,K) .GT. NUW) ,LTEST=MATCH (MATCH (I,K) -NUW, 
1IDIE2) +M+1-MOSI; GO TO 320 
LTEST=MATCH (MATCH (I,K) ,IDIM2) +M+1-M0ST 
GO TO 320 

C STORE FIRST TWO ELEMENTS OF 2-ELEMENT SET AND SET 
C POINTERS. THEN GO TO POINT WHERE NEW ELEMENT IS ADDED 
C TO SET 

365 MTEMP (1)=I 

MTEMP ii) =M ATCH/(I , K) 

INCE (1) =1 
INCH (2) =K 
GO TO 355 



C 
C 
C 
C 
C 

c 
370 
372 

375 

380 



HAVE EXHAUSTED POSSIBILITIES FOR ADDITIONS TO CURRENT SET 



OUTPUT ; 

sim a tin 

Li <1 U 



IN ANY CASE, 
LOOK FOR 



385 



405 



IS SET LENGTH OF I-NTEREST ? IF SO, 

EACH UP TO NEXT- TO- LA ST ELEMENT OF 
FURTHER POTENTIAL ADDITIONS. 

IF (M.GE.MOST) GO TO 375 
IF (M.EC.2) GO TO 405 
GO TO 380 
M 0 S T = M 

WRITE (6, 705) (MTEMP (Ml) ,MCROSS (MTEMP (Ml) ) ,M1=1,MOST) 
K = I NCR (M-1). 

J=INCR (M) + 1 
L=1 

IF (MATCH (I, K) .GT. NUW) GO 
MAT 1 = MATCH (I,K) 

MAI2=MATCH(I,IN 
L=l + 1 

IF (MATCH (MAT 1 ,L) . LT. MAT 2) 

M=M - 1 

LTEST=MATCH (MATCH (I,K) , IDIM2) + M+ 1-MOST 
GO TO 327 
K=K+ 1 



JCR (M) ) 



TO 450 



GO TO 385 



,LTEST=MATCH (MATCH (I, K) - NUW , 
0 315 



IF (MATCH (I,K) .GT. NUW) 

1IDIH2) +M+1— MOST; GO T 
LT ES1= MATCH (HATCH (I,K) ,IDIM2) +M+ 1-MOST 
GO TO 315 

C REMAINING STATEMENTS HANDLE CHECKS WHEN MATRIX 
C QUESTION IS FROM CXUW AUTO-C SUBSET 



ELEMENT IN 



415 

420 



425 



430 



MAT3=MATCH (I,K) -NUW 
L = L + 1 

IF (L.GT.LTEST) GO TO 
MAT 1 = MA1CH (M AT3 , L) 

IF (MAT 1. EQ. 0) GO TO 
MAI 2= MATCH (I ,J) -NUW 
IF (MAT2.LE.0) GO TO 
IF (MAT2..LE. MAT-1) GO 
GO TO 420 

IF (MAT I. EQ. MAT 2) , M AT 2 = M A T2 + N U W ; GO 
J = J+1 
GC TO 425 



372 

370 

370 

TO 4 30 



TO 340 
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450 MA11=MAICH (I ,K) -NUW 

MAI 2= MATCH (I , IHCR ( M) ) - NOW 

455 l= 1 + 1 

IF ^ MATCH (MAT 1,i) .LT.MAT2) GO TO 455 

LTEST= HATCH ( MA.T 1 , IDI M2 j +M+ 1-MOST 
GO TO 415 

^00 FORMAT (' 1* , ' GOOD CROSS-CORRELATICN SETS'//) 
705 FORMAT ( • O' r 5 ( IX , 7 (14 , ' ) • , 1 X,I6 ,5X) /) ) 

END 
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